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I. INTRODUCTION 

In quantum field theory symmetries and the associ- 
ated conservation laws imply Ward identities, which are 
exact relations between different types of Green func- 
tions or vertex functions^ The constraints imposed by 
Ward identities can be very useful to devise accurate ap- 
proximation schemes which do not violate conservation 
laws. In this work we shall present a non-perturbative 
derivation of the Ward identities relating the coefficients 
in the low-frequency expansion of the retarded self-energy 
E(cj + iO) of the Anderson impurity model (AIM) to 
certain thermodynamic susceptibilities. A perturbativc 
derivation of these identities has first been given by Ya- 
mada and Yosida^ For particle-hole symmetric filling 
and in the limit of an infinite bandwidth of the conduc- 
tion electron dispersion with flat density of states these 
identities read 

Re£(w + i0) = |+ (l- Xc+ 2 Xs \u + 0{u 2 ), (1.1) 



ImE(w + iO) 



. 2 2 
Xc~Xs\ u 

2A 



+ 0(w 3 ). (1.2) 



Here, U is the on-site interaction at the impurity site, 
A is the imaginary part of the hybridization function 
in the limit of an infinite bandwidth of the conduction 
electron band, and \ c and Xs are dimensionless particle 
number (charge) and spin susceptibilities, which will be 
defined in Eqs. (l4~34|) and (1435)) below. A generalization 
of Eqs. (|1.1|) and (ll.2|l to the AIM out of equilibrium can 
be found in Ref. (3J. 

Yamada and Yosida obtained the above identities by 
comparing the coefficients in the perturbation series of 
both sides of Eqs. (jl.l[) and (|1.2p to all orders in pow- 
ers of L7/A. An alternative derivation using diagram- 
matic techniques can be found in the book by Hewson. 4 
Unfortunately, in both approaches the close relation of 
the above identities to the (7(1) x (7(1) symmetry as- 
sociated with particle number and spin conservation of 



the AIM in a magnetic field is not manifest. Although 
powerful functional methods for deriving Ward identities 
non-perturbatively are well knownji apparently there ex- 
ists no derivation of Eqs. (|1.1[) and (|1 .2[) in the literature 
using these functional methods. In this work we shall 
present such a non-perturbative proof of Eqs. (|1.1[) and 
(jl.2[) by combining standard functional techniques^ with 
certain exact relations between derivatives of the self- 
energy with respect to frequency, chemical potential and 
magnetic field which we derive within the framework of 
the exact renormalization group££ 

To set the stage for our calculation and to define our 
notation, let us recall that the single-site Anderson im- 
purity model is defined in terms of the following second 
quantized Hamiltonian, 



H 



^2(e k - o-/i)4 .Cfe<7 

a 



(1.3) 



where Cka annihilates a non-interacting conduction elec- 
tron with momentum fe, energy dispersion e& and spin 
projection a, while the operator d a annihilates a local- 
ized correlated c(-electron with atomic energy Ed and spin 
projection a. The hybridization between the d-electrons 
and the conduction electrons is characterized by the hy- 
bridization energy and h is the Zeemann energy as- 
sociated with an external magnetic field. Since we are 
only interested in the correlation functions of the im- 
purity, we integrate over the conduction electrons using 
the functional integral representation of the model. 7 The 
generating functional Q c [j aj j a ] of the connected Green 
functions can then be represented as the following ratio 
of fermionic functional integrals, 



/ V[d,d]e 



-S[d„,d^]+(j,d)+(d,j) 



J V[d,d}e~ s old, 



(1.4) 
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with the Euclidean action given by 
S[d a ,d a ] = So[d a ,d a }+ Su[d a ,d a ] 
'^2GQ 1 a (iui)d UJcr d u 



+ U / dT^(T)d t (rK(r)^(r). (1.5) 
Jo 

Here, f u — jj denotes summation over fermionic 

Matsubara frequencies iui and dr denotes integration 
over imaginary time, where /3 is the inverse temperature. 
The non-interacting Green function is 



1 



iuj - £_o. a - A a (iuj) 



where 



£o,<r = E d - /i - ah 



(1.6) 



(1.7) 



is the energy of a localized d-electron with spin projec- 
tion a relative to the chemical potential fj,, and the spin- 
dependent hybridization function is given by 



iuj — £fc + /i + ah 



(1.8) 



The Fourier transform of the Grassmann fields d^(r) in 
frequency space is defined by 



J UJ 



(1.9) 



The functional Q c [j a , j a ] in Eq. (|1.4p depends on Grass- 
mann sources j a and j a and we have introduced the fol- 
lowing notation for the source terms, 

{j,d) = / ^3uad ua , {d,j)= / ^duajw (1.10) 



II. FUNCTIONAL WARD IDENTITIES 

A. U(l) Ward identities due to particle number 
and spin conservation in a magnetic field 

The Euclidean action for the correlated impurity given 
in Eq. (|1.5p is invariant under independent global U(l) 
transformations of the fields for a given spin projection. 
This symmetry implies that the generating functional 
Qc[]a:j<j] satisfies certain functional differential equa- 
tions, so-called functional Ward identities. By taking 
functional derivatives of these relations, we shall derive 
the Yamada-Yosida identities for the self-energy. Follow- 
ing Refs. PH], we perform a local gauge transformation 
on the fermion fields in the (imaginary) time domain, 

d a (r) = e-^MdUr) , 4(r) = e ia ^d' a (r), (2.1) 



where a a (T) are arbitrary real functions. The interaction 
part Sjj of our action S = Sq+Sjj is invariant under these 
transformations, so that 



S[e ia °d' a ,e- ia «d' a ] = S[d' a ,d' lT } 
! drYX{r)[d T a a {T)]d' a {r) 



if dr I dr' 
lo Jo 



E 



xA a (T-r')d' a (r') + 0(a 2 ), 



(2.2) 



where A ff (r) = J e~ lUT A a (iuj). Using the invariance of 
the functional integral representation (|1.4p of Gc[ja,ja] 
with respect to a change of the integration variables 
d„ — > d' a , d a — > d'„, and expanding to linear order in 
the gauge factors a a (r) we obtain the desired functional 
Ward identity 1,6 For our purpose it is convenient to ex- 
press the "current terms" of this Ward identity via the 
generating functional T[d a , d„] of the irreducible vertices, 
which is obtained from Gc\jcr,ja] via a functional Legen- 
dre transformation;^ 

r&.dj = (Jj d) + (d,j) - ScLw CT ] - SofaM (2.3) 

where on the right-hand side it is understood that the 
sources j a and j a should be calculated as functions of 
the field averages d a and d a by solving the equations 

SGc[ja,ja] T 6Gc%,3o] 



d„ = 



ho 



da = 



Sja 



(2.4) 



After Fourier transformation to frequency space and 
some rearrangements analogous to those in Refs. 
we obtain the functional Ward identity 



ST 



S 2 G C 



+d 



Sd 



-d, 



Sju'a'Sju'+ 

ST 



= 0, (2.5) 



where u> is an external bosonic Matsubara frequency. 
Summing both sides of this functional equation over a' 
we obtain the functional Ward identity due to the U(l) 
symmetry associated with particle number conservation, 



E 



ST T ST 



5 2 G C 



Sju'a'Sj, 



Sdu'+u, 



U)' -\-U),a' 



. (2.6) 



To obtain the analogous U(T) Ward identity associated 
with conservation of the spin projection along the axis 
of the magnetic field, we multiply Eq. (|2.5[) by a 1 before 
summing over a', which yields 

S 2 G C 



ST 



Sju'a'Sju'+u 



ST 



Sdu,' 



uj' -\-uj,a' 



= . (2.7) 
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B. SU(2) Ward identity due to spin conservation 

In the absence of a magnetic field our action (11.51) is 
invariant under arbitrary rotations in spin space. To de- 
rive the corresponding SU(2) Ward identity, we perform 
a local rotation in spin space, 



d t (r) ^ _ 



= U(r) I W 



d'dr) 



(2.8) 



where the SU(2) matrix U(t) can be written as 



U(t) = e 



«(t) 



(2.9) 



Here, cr = [a x , a v , a z ] is the vector of Pauli matrices 
and ol(t) is a time-dependent three-component vector. 
Expanding to linear order in a(r) we obtain, after the 
same manipulations as in Sec. Ill Al the following SU (2) 
Ward identity, 



s 2 g c 



^Juj'(T / ^juj'-\-UJ.G 



ST 



^<^o/-|-cj,ct 



ST 



where the superscript i = x,y,z labels the three compo- 
nents of the vector operator cr. Together with the particle 
number conservation Ward identity (|2.6j) these equations 
are equivalent to the four Ward identities 



5T 



6 2 G C 



5T 



0, (2.11) 



where a, a' € {t, -1} are now fixed spin projections. Note 
that in the absence of a magnetic field the Green function 
and the self-energy are independent of the spin quantum 
number cr, so that we may write G a (iuj) — G(iu) and 
Y, a {iuj) — S(iw). For i — z we have <j% a , = o'8 a .<T' so 
that Eq. (|2.10[) reduces to the zero-field limit of the U(\) 
Ward identity (I2.7[) associated with the conservation of 
the spin component in the direction of the magnetic field. 



III. WARD IDENTITIES FOR THE 
SELF-ENERGY 

A. Particle number conservation 

For the AIM in a magnetic field, the first terms in the 
functional Taylor expansion of the generating functional 



T[d a , d a ] are of the form 6 

rfJcr,^] = r + / Efffa)^^ 

Jul u 



- ''-"1 era' 



+0 (cFc^ 



(3.1) 



where To is proportional to the interaction correction to 
the grand canonical potential, £ CT (iw) is the exact irre- 
ducible self-energy, and Uy^, {iu}' Xl iuj' 2 \ iui2, iu>i) is the ex- 
act effective interaction. The term in the second line of 
our functional Ward identity (12.61) then yields 



E 



din 1 n' ' 



8T 



duj' 



ST 



r- j ^LU' -\-OJ.(j' ' r- j 

OCLuj' ,(7 / O&tjji a > 

u 

= / y^}T, a '(iui' + iui) - Eg^iui'^daj'+ua'du,'^ 

+0 (cPd 2 ) . (3.2) 

To extract the Ward identity for the self-energy from our 
functional Ward identity (12.6[) . we take the second func- 



0, (2-10) tional derivative 



of both sides of Eq. 



and set then all fields equal to zero, which yields 

T, a (ioj + iuj) — Y, a (iuj) 

£fe(«j' + ifl)-Go>(iw') 

>' a' 

xG„>(iuj' + iu))G a >(iu:') 
xC/^o-' (* w + + 



(3.3) 



To obtain the right-hand side of this Ward identity, we 
have used 



s 2 g c 



ficlds=0 



tidier Sd-jj-^/jj^ (^Jcj'CT'^Joj'+a 

G al {iJ)G al {iuj' + iuj) 

x V^\, (ilu + iu), iuj'; iuj 1 + iuj, iu), 



(3.4) 



which follows from the tree expansion relating the con- 
nected Green functions generated by Q c to the irreducible 
vertices. In the limit of vanishing bosonic frequency 
Co — > the Ward identity (|3 . 3[) reduces to 



E 



1 



dA a ,(lLU>) 



d{iuj') 



(3.5) 



d(iu>) 

xT a ^(iuj,iuj'), 

where we have defined 

T a ^ <jl (iuj 1 iuj') = Uyl,(iu>,w']iu>\kj). (3.6) 

The above identities are valid for arbitrary hybridiza- 
tion functions A a (iuj). Of special interest is the limit of 
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infinite bandwidth of the conduction electron band with 
flat density of states, where the general expression for 
A a (iu}) given in Eq. (jl.8|) is given by 



A a (iuj) = -iAsgnw = -iA[26(w) - 1]. 



(3.7) 



Here, the hybridization energy A is assumed to be inde- 
pendent of the chemical potential and the magnetic field. 
In this limit 



dA a ,{iJ) 
d(iui') 



= -2A5(w'). 



(3.8) 



The term (dA a ,(iuj')/d{iu'))Gl,(iu') in Eq. (J33J) is then 
ambiguous because the (5-function in the first term is 
multiplied by the sign-function associated with the hy- 
bridization function in G^,(iui'). To properly define this 
term one should use the Morris-Lemmaj^Jii which states 
that the product of the delta-function with an arbitrary 
function f(Q(x)) of the 0-function should be defined via 



5(x)f(Q(x)) = S(x) / dtf(t). (3.9) 
Jo 

We conclude that in the infinite bandwidth limit 
- d -^0yGl^) = 2A5{J)Gl,{ l J) 

r 1 i 

= 2A<J(w') / dt— -— -pr 

V 'Jo K ff '+«A(2t-l)] 2 

2A 

= <V)^-j^ = 2kS("')pA0), (3-10) 

where £ CT = £o,cr + S CT (iO) is the true excitation energy of 
a localized electron with spin projection a, and 



A 1 

MO) = — c2 i A2 
7T & + A 2 



(3.11) 



is the spectral density of the localized d-electron with 
spin a at vanishing energy. Substituting Eq. (|3T0|) into 
p.5[) we obtain the corresponding Ward identity in the 
infinite bandwidth limit ^ 

-^^(Ojr^Ciw.O). (3.12) 



If we (incorrectly) replace Pcr'(0) - ► Po-(O) in the second 
line of Eq. (|3T2]l . we arrive at Eq. (5.71) of Ref . 0]. 



corresponding functional Ward identity (|2.7I) . By com- 
paring the functional Ward identity (|2.6j) due to particle 
number conservation with the corresponding U(l) Ward 
identity (|2.7p due to spin conservation, we conclude that 
in the spin case the Ward identities for the self-energy 
can be obtained from those of the particle number case 
by simply replacing 



(3.13) 
(3.14) 



In particular, the spin-analogue of the particle number 
Ward identity ([3~5]l is 



d(ioj) 



i 



dA a ,{iu') 



xTo-.o-' (iw, iuj'). 



G\> (iu)') 
(3.15) 



Combining this equation with the corresponding equa- 
tion for the particle number, Eq. p.5[) , we obtain the 
two Ward identities 



9E CT (iw) 
d(iuj) 



1 



dA a (iw') 



d{iLu') 
dA_ a (iu' 



(3.16) 



d{iuj') 



G 2 _ CT K) 



(3.17) 



In the limit of an infinite bandwidth these equations can 
be written in analogy to Eq. (|3.12p as 



d(iui) 



Gl (iw')r CTiff (io;, io/) 

M^o)r CT ^(iw,o), 



(3.18) 

-p_ CT (0)r CT ,_ CT (ia;,0). (3.19) 

Finally, we note that in the absence of a magnetic field 
the SU (2) functional Ward identity (|2.10[) does not imply 
any further independent relation for the self-energy, since 
in this case we just have H a (iw) = £(«<;). However, 
taking higher order functional derivatives of Eq. (|2.10j) 
we may obtain symmetry relations involving higher order 
vertices, which are beyond the scope of this work. 



IV. PROOF OF THE YAMADA-YOSIDA 
IDENTITIES 



B. Spin conservation 



A. Relations between derivatives of the self-energy 



The Ward identity for the self-energy associated with 
the conservation of the spin projection in the direction of 
the magnetic field can be derived analogously from the 



The above Ward identities can now be used to derive 
exact relations between the derivatives of the self-energy 
with respect to frequency, the chemical potential and the 



5 



magnetic field. In the book by Hewson^ one can find a 
derivation of these identities using diagrammatic argu- 
ments. Here, we show that the correct relations can be 
obtained quite simply within the framework of the ex- 
act renormalization group (also known as the functional 
renormalization group). 5,6 As a special case of the gen- 
eral exact renormalization group equation for the irre- 
ducible self- e nerg y of an interacting Fermi system derived 
in Refs. [6|fT^fl3| we find for the self-energy of the AIM 



<9a£<x(«w) 



(4.1) 



where we have used again the notation (|3.6p for the ef- 
fective interaction, and where A is some flow parameter 
(cutoff) appearing in the Gaussian part of the action. 
All Green functions and vertices in Eq. (|4.1[) implicitly 
depend on the parameter A. The so-called single-scale 
propagator is defined by 



G a {vJ) = -Gl(iu:)d A Gu l(iuj). 



(4.2) 



The exact renormalization group flow equation (14.11) is 
valid for any choice of the flow parameter A. In partic- 
ular, we may choose the chemical potential fi as a flow 
parameter^ Then A = fi and hence 



d 

dA G o,a( iuj ) = -g^[ iu} - E d + M + ah 



A ff (iw)] 



1 



dA a {uj) 
dfi 



so that the single-scale propagator is simply 

dA a {iuj) 



G CT M = -G£(iw) 



1 - 



(4.3) 



(4.4) 



In this chemical potential cutoff scheme the exact renor- 
malization group flow equation (|4.1[) reduces to 



dfi 



E 



i - 



Of i 



Gl,(iu') 



(4.5) 



To further manipulate this expression, let us note that for 
u) 7^ the integrand of the hybridization function given 
in Eq. (|1.8[) is non-singular, such that 



dA a (iuj) dA„(iu>) 



d(ioj) 



dfi 



LU ± 0. 



(4.6) 



For w = 0, however, we have to be more careful. Defining 
the spectral density 



fJ 



(e)=nJ2\V k \ 2 6(e-e k ), 



(4.7) 



Using the well-known formula 



1 1 . 

x ± iQ + x 



(4.9) 



where P is the principle value, we see that A a (z) has a 
branch cut along the real axis. More explicitly, defining 



A = g(n + ah), 



we obtain 



dA a {iuj) dA a (iui) 



d(iuj) 



2AS(uj), 



(4.10) 



(4.11) 



where the term — 2A5(uj) arises from the discontinuity 
across the branch cut. In the limit of a flat band of 
infinite width discussed above, g(e) = A for all e and 
dA a (iuj)/dfi = 0, so that Eq. (|4"TTj) reduces to Eq. (|3~5|) . 
Using the identity (|4.11[) we may write our exact renor- 
malization group flow equation (|4.5[) in the form 



dfi 



E 



d{iu') 



xGl,(icj')T a ^(iu},iu)') 



E 



dA a ,(iu') 



d(iuj r 

^p CT <(o)r CT ^( iW ,o) 



Gl,(iLo')T a ^,(iuj, faj') 



(4.12) 



Comparing the right-hand side of this exact relation with 
the right-hand side of the particle number Ward identity 
(|3.5[) we conclude that for all frequencies we have the 
following exact identity, 



dT, a (iLo) dT, a (iuj) 
d(iuj) dfi 



$^p«r/(0)IV(iu;,0). (4.13) 



Next, let us choose in our exact renormalization group 
flow equation (|4.ip the magnetic field as a flow parameter 
(A = h). The relevant single-scale propagator is then 



G a (itu) = -Gl(iu) 



dA^iu) 



and hence 

dT, a {iui) 
dh 



-- X' 



1 - a 



Oh 



,0Av(iw') 



dh 



xT a<(T >(iui, ioj'). 



(4.14) 



Gl,{iJ) 
(4.15) 



the hybridization function A a {iuj) can be rewritten as 

1 



1 f°° 
A a (ioj) = — de 

K J-oo 



9(e) 



e + jU + ah 



(4.8) 



Noting that 



_A a (iuj) A a (ioj) 



dh 



(4.16) 
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Eq. (|4.15l) can also be written as 



E*' 



<9(iu/) 

^£rVa'(0)r CT) a'(^.Q) 



G^(iw')r o - !CT '(iw, iu/) 



(4.17) 



Comparing this with the f7(l) spin Ward identity (|3.15|1 . 
we conclude that 



dH a (ioj) dTi a (iuj) 
= a- 



d(iu) 



dh 



- a ^ o'p*' (0)r a , a /(iw,0). 

(4.18) 



Of particular interest are the above relations for u) — > 
0, because in this limit the w-derivative determines the 
wave-function renormalization factor Z„ via 



d{iu) 



i-z-\ 



(4.19) 



j=0 



Taking the limit u> ->• in Eqs. (|4~T3| and (|4~18| and 
using the fact that due to antisymmetry the effective in- 
teraction at vanishing frequencies has the form 

U ( a % (0, 0; 0, 0) = IV, (0, 0) = 5 CT ,_a'IV CTJ (4.20) 
we obtain 



d{iuj) 



d(iuj) 



j=0 



j=0 



dp 

as ff (z0) 



- p_ a (p)T„ t - a , (4.21) 
+ ^(0)1^.(4.22) 



Adding and subtracting these equations we obtain 



d(iu) 

2p_ ff (o)r 



1 cr — , (4.23) 



j=0 



9/i <9/i 
9S CT (iO) 9S CT (i0) 



(4.24) 



9/i i% 
Eq. (|4~2li|) is a corrected version of Eq. (5.80) of Ref.0]. 

B. Yamada-Yosida identities 

To make contact with the work of Yamada and 
Yosida£, we now assume a flat density of states and take 
the limit of an infinite width of the conduction electron 
band, D — > oo. Then there is a simple exact relation be- 
tween the occupation numbers n a of the impurity level 
and the self-energies T, a (i0) at vanishing frequency, 4 



1 1 

n„ = arctan 

2 7T 



Ed — p — oh + S ff (i0) 



(4.25) 



Taking derivatives of this expression with respect to p 
and h we obtain 



W = ^ (0) 



1 - 



1 -CT 



dp 

as a (z0) 



a/?. 



Hence, 



as g (io) 

dp 
dh 



= 1 - 



= 1 - 



1 drier 
p CT (0) <9/z ' 

1 dn a 
p a (0) a ~dh 



(4.26) 
(4.27) 

(4.28) 
(4.29) 



Substituting these relations into the Ward identities 
(|4.23[) and (|4.24[) we obtain for the w-derivative of the 
self-energy 



d(iuj) 



1 



1 



2p CT (0) 



dp 



_ dn a 
dh 



and for the effective interaction at zero energy 



i 



2pt(0)pi(0) 



dn a 
dp 



_ dn a 
~dh 



(4.30) 



(4.31) 



For h — > the self-energy and the density of states are 
independent of the spin projection, so that we may write 
E CT (iw) = E(tw), Pa{0) = p(0), and r CT ,_ CT = F ± . As 
in Ref. 0, we now focus on the particle-hole symmet- 
ric case, where p(0) = 1/(ttA) in the flat band infinite 
bandwidth limit considered here. Averaging both sides 
of Eqs. (|4.30|) and (|4.31[) over both spin projections, we 
obtain the Ward identities 



d{ioj) 



p(o)rj 



Xc + Xs 

2 



(4.32) 
(4.33) 



where the dimensionless particle number (charge) and 
spin susceptibilities are defined by 



Xc 



_ Xc _ \ "» VUg 



p(o) 



7rA 

T 

7rA 



9? 



dn a 
~dh 



(4.34) 
(4.35) 



The identities (I4.32[) and (|4.33[) have first been obtained 
by Yamada and Yosida^ who showed that both sides of 
these equations have identical series expansions in pow- 
ers of f//A. Such a perturbative proof relies on the as- 
sumption that there are no non-analytic terms which are 
missed by the series expansions. Our proof of Eqs. (|4.32[) 
and (|4.33[) given above shows more clearly that these 
identities are a direct consequence of the U(l) x U(l) 
symmetry associated with particle number and spin con- 
servation of the AIM in a magnetic field. Moreover, our 
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derivation of Eqs. (|4.32[) and (|4.33l) is non-perturbative, 
because it relies only on the symmetries of the AIM and 
the associated functional Ward identities, and on the 
exact renormalization group flow equation for the irre- 
ducible self-energy. 

Using relation ()4.19|) . we see that the identity (|4.32p 
implies that the wave-function renormalization factor of 
the AIM can be expressed in terms of the susceptibilities 



(4.36) 



Xc + Xs 



The other Ward identity (|4.33j) allows us to express the 
imaginary part of the retarded self-energy to the suscep- 
tibilities. Therefore we recall that the skeleton equation 
for the self-energy of the AIM (which is a consequence 
of the Dyson- Schwinger equation, which in turn follows 
within a functional integral approach from the invariance 
of the functional integral under infinitesimal shifts of the 
fields^) implies the following exact expression for the 
imaginary part of the self-energy of the symmetric AIM 



ImE(w + iO) 



(4.37) 



Substituting the identity (|4.33[) for the dimensionless ef- 
fective interaction F_l/(7tA) = p(0)T±, we obtain 



ImS(a; + iO) 



2 2 

Xc ~Xs\ 

2A 



+ 0(u; 3 ). (4.38) 



Taking the different notations for the dimensionless sus- 
ceptibilities into account^ Eq. (|4.39p agrees with the 
expressions derived by Yamada and Yosida.— 



V. CONCLUSIONS 

In this work we have used modern functional meth- 
ods to give a non-perturbative proof of the Yamada- 
Yosida identities, which express the coefficients in the 
low-frequency expansion of the self-energy of the An- 
derson impurity model in terms of thermodynamic sus- 
ceptibilities. In contrast to the derivation of these rela- 
tions given by Yamada and Yosida^ which is based on 
a series expansion in powers of the interaction, our non- 
perturbative proof relies on exact Ward identities and on 
an exact renormalization group flow equation for the irre- 
ducible self-energy. From our derivation it is obvious that 
the Yamada- Yosida identities are a direct consequence 
of the U{\) x U(l) symmetry of the Anderson impurity 
model associated with the conservation of the particle 
number and the total spin component in the direction of 
an external magnetic field. 

We have also presented a thorough derivation of the 
general functional Ward identities of the AIM due to 
particle number and spin conservation. Furthermore, we 
have shown that various identities relating the derivatives 
of the self-energy with respect to frequency, chemical po- 
tential, and magnetic field can be obtained by combin- 
ing the Ward identities with exact renormalization group 
flow equations for the self-energy. 



On the imaginary frequency axis, the low-frequency ex- 
pansion of the self-energy of the symmetric AIM in the 
infinite bandwidth limit is therefore 



SH = | + ( ' 



Xc+Xs\ . , .{Xc-Xs) 2 2 



iuj + i- 



-uj sgnu; 



2 j— ' ■ 8 a 
analytic terms 0(uj 2 ). (4.39) 
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